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J3 ■ Abstract 

We investigate the theoretical uncertainties of P-wave charmonium decays 
^ ■ into two pions, Xcj ~^ ir°n . Constraining the pion distribution- 

amplitude from the recent precise data on F nl (Q 2 ), we find an order-of- 
magnitude discrepancy between data and prediction. The disagreement 
persists even after inclusion of transverse degrees of freedom and Sudakov 
suppressions. We propose the colour-octet mechanism as the solution to the 
puzzle. The colour-octet decay contribution arising from the higher Fock 
component |ccg) of the Xcj wave function is actually not power suppressed 
with respect to the usual colour-singlet decay arising from the dominant 
|cc) Fock state. An explicit calculation yields an agreement with the data 
for a very reasonable value for the single extra non-perturbative parameter. 
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Up to very recently, hard exclusive reactions involving pions could be believed to be 
successfully described within the hard-scattering approach. The reactions to be mentioned 
are the pion-photon transition form factor, the pion electromagnetic form factor, pion-pair 
production by two photons, as well as the two-pion decays of the charmonium states J/ip, 
Xco, and Xc2- The hard-scattering approach (HSA) provides the framework of calculations 
of exclusive reactions at large momentum transfers a full amplitude factors into 
two parts, a hard-scattering amplitude T#, calculable in perturbative QCD, and parton- 
distribution amplitudes <3># for each hadron H. The amplitude T# describes the scattering 
of clusters of collinear partons from the hadron. The leading-twist contribution, i.e. the 
one with the weakest fall-off with Q, the typical momentum transferred in the process, is 
given by valence-parton scatterings only. Hence the only non-perturbative input required 
are the probability amplitudes (or quark distribution amplitudes) $H(xi,Q) for finding 
valence quarks in the hadron, each carrying some fraction Xi of the hadron's momentum. 

In the case of mesons, the leading Fock state is the |qq ) state, which, for the pion in the 
isotopic limit, can be described by a symmetric wave function $> w (x) = $^(1 — x). Upon 
expansion over Gegenbauer polynomials C^l 2 \ the distribution amplitude at a scale \if 
can be characterized by non-perturbative coefficients B n (see and references therein): 



1+ £ bJ^Y C^\2x-l) 

n=2,4,... \<W0)/ 



Here f n is the pion decay constant, f n = 130.7 MeV, a s the strong coupling constant, and 
/io a hadronic scale, 0.5 £ fio £ 1 GeV. Since the 7„ in (P are positive fractional num- 
bers increasing with n, higher-order terms are gradually suppressed and any distribution 
amplitude evolves into (pas(x) = 6x(l — x) asymptotically, i.e. for ln(Q 2 /Aq CD ) — > oo. 

Lowest-order perturbative QCD (pQCD) calculations in leading-twist of the above- 
mentioned reactions involving pions were all in reasonable agreement with data, provided 
the quark distribution amplitude in the pion is strongly end-point-concentrated (in sharp 
contrast to the asymptotic distribution amplitude). Such a distribution amplitude is given 
by the Chernyak-Zhitnitsky (CZ) one defined by B 2 = 2/3 and B n = for n > 2. 
So far, there existed basically only one exception, namely the pion form factor in the 
time-like region: the HSA prediction with the CZ distribution amplitude fails by about 
a factor of 2 as compared with the data. The experimental information on the time-like 
form factor comes from two sources, e + e~ — > 7r + 7r~ and J/ip — > 7r + 7r _ , which provides, to 
a very good approximation, the form factor at s — M],^. Although the e + e~ annihilation 
data of Ref. || suffer from low statistics, they agree very well with the result obtained 
from the J/ip decay. 

Meanwhile, the situation has changed drastically. On the experimental side, very pre- 
cise new data for fairly large momentum transfer || have shown that F 1T1 {Q 2 ) is smaller 
than predicted from the CZ distribution amplitude. On the theoretical side, the HSA has 
been modified through the incorporation of transverse degrees of freedom and Sudakov 
suppressions . This more refined treatment (termed the "modified HSA" ) allows the 
calculation of the genuinely perturbative contribution self-consistently, in the sense that 
the bulk of the perturbative contribution is accumulated in regions of reasonably small 
values of the strong coupling constant a s . This approach hence overcomes arguments 0] 
against the applicability of the standard HSA (sHSA: the lowest-order pQCD approach in 
the collinear approximation using valence Fock states only) to experimentally accessible 
regions of momentum transfers. It is to be stressed that the effects of the transverse 
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degrees of freedom taken into account in the modified HSA (mHSA) represent soft con- 
tributions of higher-twist type. Still, mHSA calculations are restricted to the dominant 
(valence) Fock state. 

A detailed analysis of the photon-pion transition form factor F 1X1 (Q 2 ) within the 
mHSA H [| shows that the pion distribution amplitude is close to the asymptotic form 
(f)as(x). In order to give a quantitative estimate of the allowed deviations from the asymp- 
totic distribution amplitude, one may assume that B 2 is the only non-zero expansion 
coefficient in ([!]). The truncated series suffices to parametrize small deviations. More- 
over, it has the advantage of interpolating smoothly between the asymptotic distribution 
amplitude and the frequently used CZ distribution amplitude. For large momentum 
transfer our assumption is justified by the properties of the anomalous dimensions 7„. 
Quantitatively, choosing /i = 0.5 GeV, a best fit to the F nl data above 1 GeV 2 yields 

B 2 {fio) = -0.006 ± 0.014 mHSA 

= -0.39 ± 0.05 sHSA . (2) 

The non-perturbative pion distribution amplitude (i.e. the coefficient B 2 ) is, in fact, 
best extracted from F 7r7 (Q 2 ), since the pseudoscalar meson-photon transition form fac- 
tors are pure QED processes with QCD corrections of only 10-20% and higher Fock-state 
contributions suppressed by powers of a s /Q 2 and expected to be small too |jI0| . Uncer- 



tainties due to these effects may be absorbed into a theoretical error of B 2 , which we 
generously estimate to amount to ±0.1 for the mHSA and ±0.39 for the sHSA. However, 
with a distribution amplitude close to the asymptotic form, calculations for the other 
processes involving pions disagree strongly with data, at least when calculated in the 
collinear approximation. In this paper we shall investigate in detail the Xcj decays into 
two pions, by calculating the decay rates and branching ratios in the modified HSA, using 
the information on the process-independent expansion coefficient B 2 obtained from the 
7T7 form factor. 



Using the method proposed by Duncan and Mueller |TT[, the decay rate calculated in 
the sHSA can be written as 

T[ Xc j^n + n-] = f^\R' P (0)\ 2 a s (fiT 1 ) 4 a J + bj B 2 (fi ) + cj B 2 (fi ) 2 2 , (3) 

where aj, bj, and cj are (analytically) calculable real numbers, see table [TJ. The evolution 
factor [ln(/i 2 7 /AQ CD )/ln(/iQ/AQ CD )] 72 (see ([[])) is absorbed into the coefficients bj and cj 
(quadratically in the latter case), and 72 = 50/81. //ff d is an appropriate renormalization 
scale, which is of the order of the charm-quark mass m c [I], |I2, [H| . In the sHSA analysis, 



it is customary to identify factorization scale hf and renormalization scale. Finally, R'p(0) 
denotes the derivative of the non-relativistic cc wave function at the origin (in coordinate 
space) appropriate for the dominant Fock state of the Xcj, a cc pair in a colour-singlet 
state with quantum numbers 2S+1 Lj = 3 Pj. 

In the following we will choose as central values -Rp(O) = 0.22 GeV 5//2 and m c = 
1.5 GeV, which is consistent with a global fit of charmonium parameters |T4| as well as 



results for charmonium radii from potential models[T5|]. With this choice, /x^ = m c , 
and B 2 (fi ) = —0.39 (see (^)) the sHSA results for the Xc0{2) —> 7r + 7r~ decay widths are 
almost two orders of magnitude below the experimental data, see table |[ To assess the 
uncertainty of the result we vary the charm-quark mass between 1.35 GeV and 1.8 GeV. 
Of course, the value of R' P (0) has to be adjusted accordingly: using the well-known 



2 



Approach 


a 


b 


c 


sHSA 


19.60 


37.43 


26.23 


mHSA 


-14.67 + 23.36 % 


15.36 + 73.22^ 


27.18 + 26.39 2 




a 2 




c 2 


sHSA 


4.589 


8.479 


7.719 


mHSA 


-3.420 + 5.1412 


6.681 + 13.42 2 


5.744 + 1.3842 



Table 1: Coefficients defined in obtained within the standard HSA (sHSA) for //ff d = 
m c = 1.5 GeV and coefficients ( |i~0l) of the modified HSA (mHSA). 



B 2 


m c 


Aqcd 


T(XcJ -> vr+TT-) [keV] 


BR( Xc J -> 7T + 7T-) [%] 




[GeV] 


[GeV] 


J = 


J = 2 


J = 


J = 2 


Standard HSA 


-0.39 


1.5 


0.2 


0.872 


0.065 


0.011 


0.003 





1.35 


0.25 


15.3 


0.841 


0.113 


0.020 


Modified HSA 





1.5 


0.2 


8.22 


0.41 


0.102 


0.017 


0.1 


1.5 


0.2 


12.13 


0.53 


0.151 


0.021 


-0.1 


1.5 


0.2 


5.61 


0.33 


0.070 


0.013 





1.8 


0.2 


2.54 


0.12 


0.050 


0.008 





1.35 


0.2 


15.3 


0.78 


0.144 


0.024 





1.5 


0.15 


4.34 


0.21 


0.071 


0.011 





1.5 


0.25 


13.1 


0.68 


0.128 


0.022 


Standard HSA plus colour-octet contributions 





1.5 


0.2 


49.85 


3.54 


0.36 


0.177 


Experiment 


PDG[16 




105 ± 30 


3.8 + 2.0 


0.75 + 0.21 


0.19 + 0.10 


BES[17| 




62.3 + 17.3 


3.04 + 0.73 


0.427 + 0.064 


0.152 + 0.034 



Table 2: Decay widths and branching ratios of XcJ — * vr + 7r~ for various choices of the 
parameters compared with experimental data. The standard HSA results are obtained 
with //|f = m c , the branching ratios are evaluated with = m c and F± = 1. 
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Figure 1: Feynman diagrams for the colour-singlet decay XcJ — » 7r 7r ( J = 0, 2). 

scaling properties of quarkonium potential models [IS] we take R' P (0) = 0.274 GeV 5//2 
for m c = 1.8 GeV and R' p (0) = 0.194 GeV 5 / 2 for m c = 1.35 GeV. However, even when 
stretching all parameters to their extreme values, the predictions stay a factor 3-6 below 
the data. Note thatfl a s (m c ) = 0.447 for m c = 1.35 GeV and A QCD = 0.25 GeV. 

The starting point of the calculation of the Xcj decays within the modified HSA is the 
convolution formula 



32 y/2 7r 3 / 2 R ' p (0) 



M( X cj -> vr 7T ) = : a j 

irr 

, r d 2 b 



Jo dxdy J 7^)2 ^^' b ) Tf/J ( x ' 2/ ' b )* 7r ( x ' b ) ex v[-S(x,y,b,t u t 2 )} , (4) 

which adapts the methods proposed by Botts and Sterman to our case (o~o = 1; 
a 2 = v^3/2 ). Note that we work in the transverse coordinate space (the quark-ant iquark 
separation b is canonically conjugated to the usual transverse momentum k±). Thj is the 
Fourier-transformed hard-scattering amplitude to be calculated from the graphs shown in 
Fig. ft 

T HJ {x, y, k ± -k ± ) = + ^ ^ + k) N ( 1 + — — 1 , (5) 



and with k = (kj_ — k'j_)/2m c 



iV = s(l - j/) + (1 - + 2 k 2 , 

2 i 2 

g 1 =xy-k , 

fl g = (l- a .)(l-y)-k a . (6) 



In analogy to previous applications of the modified perturbative approach, the Sudakov 
exponent is given by 

S(x, y, b, tx, t 2 ) = s(x, b, 2m c ) + s(l - x, b, 2m c ) + s(y, b, 2m c ) + s(l - y, b, 2m c ) 
4 log(t 1 /A QCD )log(t 2 /A Q c D ) 

/3 ° g log(l/(6A QCD ))2 ' {<) 



1 We calculate a s in the one- loop approximation with Aqcd = 200 MeV and n/ = 4. 



where the function s(x, 6, Q) can be found, for instance, in Ref. fT9|| . 

The renormalization scales tj = max(2 v /x^'m c , 1/6) (i = 1, 2) appearing in a s and in 
the Sudakov exponent are determined by the virtualities of the intermediate gluons, which 
depend non-trivially on the integration variables. This choice of the renormalization scale 
avoids large logs from higher-order pQCD. The factorization scale is given by the quark- 
antiquark separation 6, fip = 1/6 = 1/|6|. The ratio 1/6 marks the interface between 
non-perturbatively soft momenta, which are implicitly accounted for in the pion wave 
function and the contributions from semi-hard gluons, incorporated in a perturbative 
way in the Sudakov factor. 

The last object appearing in (4) is the full (soft) wave function of the pion describing 



also the dependence on the transverse momentum. Following |20| we make the ansatz 



^ n (x,b;fx F ) 



2^6 
^(x, 6) = An exp 



x(l — x)b 2 



4a2 



Note, however, that a n is not a free parameter since it is fixed from tt° — > 77 | ZIfl . That 
constraint leads to the closed formula l/a% = 8 (1 + B 2 ) n 2 f 2 . This additional B 2 
dependence is taken into account in our calculation by expanding the Gaussian in E over 
B 2 . For the asymptotic form of the wave function (B 2 = 0) the transverse size parameter 
a n takes the value 0.861 GeV -1 . 

In terms of the amplitude (4) the decay widths are given by 



r(x c o 
r(x C 2 



71 71 



7T 7T 



32 n m c 
1 

240 7r m c 



M( X co - 

M{ X c2 



71 71 



(9) 



The final results in the modified HSA, tables [1] and ^] can be obtained numerically only, 
but can be cast into a form similar to (|3T) 



nXcJ^n+n-] = f* — \R' P (0){' 



77^ 



a j + 6 j B 2 {ix ) + cj B 2 ([i C 



(10) 



The coefficients aj, bj, and cj are now complex- valued. It is still convenient to divide 
out the fourth power of a s at the fixed scale m c in (|10D , since the main effect of the 
strong coupling is thus made explicit. The actual effective renormalization scale fi p in 
the modified HSA differs from m c and depends on B 2 . We find 



/4 ff 



1.15 GeV 
0.85 GeV 



off 



0.43 (B 2 = 0) 
0.52 (B 2 = 2/3) . 



Note that part of the increase of the decay widths compared to the sHSA results follows 
from the larger values of the strong coupling constant a s , cf. ( |Tl"D with a s {m 2 ) = 0.374. 

The calculation of the charmonium decay is theoretically self-consistent in so far as 
only 2% (20%) for B 2 = (B 2 = 2/3) come from soft regions where the use of pQCD is 
inconsistent. The soft region is defined by a s (t1) a s (t 2 ) > 0.5. 
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Uncertainties in our predictions of the Xc decay rates into pions arise from uncertainties 
in the values of B 2 , Aqcd> -Rp(O) and the charm-quark mass m c . In order to show the 
dependence on m c we give our results in table |2] for the above three choices of m c (with 
R'p(0) adjusted accordingly). We also show the dependences on Aq CD and B 2 , where we 
stretch B 2 up to the maximal value allowed by the analysis of the photon-pion transition 
form factor, see above. As shown in table || the predicted decay widths for both Xco an d 
Xc2 are well below the data, even if the parameters are pushed to their extreme values. 

In order to eliminate the dependence on R'p(0) and reduce the one on m c one can con- 
sider the branching ratios into tt + tt~ rather than the absolute widths. When normalising 
to the total width we correct for the fraction pj of Xcj decays that goes into light hadrons 
(rather than decaying electromagnetically or into lighter charmonium states, and which 
can be determined from data). The branching ratio has a form similar to fllPl): 



BR[ XcJ -> 7T+7T-] = F 1 7 Jj ^ dj aj + bj B 2 (im>) + cj B 2 (p 

a s {n % pY mi 



2 



(12) 



where d = po/6, d 2 = 5p 2 /8, and Fx = 1 (see below). The results for the branching 
ratios, see table 0, are indeed less sensitive to the parameter choice (and independent of 
R' P (0)), and confirm our conclusion based on the decay widths: The calculation based on 
the assumption that the Xcj is a pure cc state, is not sufficient to explain the observed 
rates. The necessary corrections would have to be larger than the leading terms. A new 
mechanism is therefore called for. 

In the quark-potential model, any charmonium state is a pure cc state, where the 
quarks have angular momentum L and spin S such that these match the quantum numbers 
by which any meson is characterized in QCD, total spin J, parity P and charge conjugation 
C: J = L + S, P = (-1) L+1 , C = (-1) L+S . The Xcj states are therefore made out of 3 Pj 
cc states (in the spectroscopic notation ( 2S+1 ^Lj), which, obviously, are in colour-singlet 
(c = 1) states. Although the HSA acknowledges higher Fock components in the meson's 
wave function, these do not contribute in either sHSA or mHSA calculations since their 
contributions are assumed to be suppressed by powers of the hard scale JTJ. 

Recently, the importance of higher Fock states in understanding the production and 



the inclusive decays of charmonia has been pointed out [I3|. The heavy-quark mass 
allows for a systematic expansion of both the quarkonium state and the hard, short- 
distance process and, hence, of the inclusive decay rate or the production cross section. 
The expansion parameter is provided by the velocity v of the heavy quark inside the 
meson. In the case of the Xcj, the Fock-state expansion starts as 

|Xcj) = 0(1) \ccx( 3 Pj)) + 0(v) \cc 8 ( 3 Sx) g) + 0(v 2 ) (13) 



where the subscript c specifies whether the cc c is in a colour-singlet (c = 1) or colour-octet 
(c = 8) state. 

The crucial observation is now that, for inclusive decays (and also production rates), 
both states in (^) contribute at the same order in v and, hence, the inclusion of the 
"octet mechanism" , i.e. the contribution from the |ccgg) state, is necessary for a consistent 
description. (Without its inclusion the factorization of the decay width into long- and 
short-distance factors is spoiled by the presence of infrared-sensitive logarithms.) In simple 
terms this can be understood by realizing that the annihilation of the cc pair of the higher 
Fock state into light hadrons is an S'-wave process compared to the P-wave annihilation 
of the leading Fock state. The inclusion of the colour-octet contribution to the total Xcj 
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Group 3 




Group 5 




Group 6 




For decay into 
neutral pions 



Group 11 



Figure 2: Representatives of the various groups of colour-octet decay graphs. 



width does, however, worsen the comparison of flT2|) with data since the fraction Fi of the 
total hadronic width, which proceeds through the c = 1 contribution, is less than unity. 

Here we propose the inclusion of contributions from the |ccg( 3 S'i)g) Fock state to the 
exclusive Xcj decays as the solution to the failure of the HSA. We shall show that this 
colour-octet contribution is indeed not suppressed by powers of l/m c compared with the 
conventional, colour-singlet contribution. Let us first, however, discuss the problem of 
colour conservation. This is a new feature, specific to exclusive decays. 

The obvious solution to colour conservation seems to be to demand one of the final- 
state pions to be also in a higher Fock state, i.e. to consider |ccg) — ► |qqg) + |qq). The 
hard process consists of the Feynman diagrams of Fig. [I] plus the diagram in which c and 
c annihilate into a single gluon, which, in turn, makes a qq pair. Yet, we disregard this 
possibility. The reason has to do with the importance of higher Fock states in the pion. 
Since the qq pair of the dominant Fock state is in a 1 S state, it is in a 1 P\ state for the 
next-higher state. The first S*-wave state |qq~8( 3 >S'i)g) is reached at 0(v 2 ). Contributions 
from higher Fock states of the pion therefore seem to be suppressed. 

Our solution to colour conservation is the following: We attach the gluon of the |ccgg) 
state (in all possible ways) to the hard process leading to the Feynman diagrams shown in 
Fig. |2|. We do, however, realize that the virtualities of the quark propagators to which the 
soft gluon couples are typically (much) smaller than the virtualities associated with the 
hard gluons. Therefore we treat the coupling of the soft gluon as a parameter separate 
from a s and denote it by a s s oit . It will turn out that the final result depends on just a 
single non-perturbative parameter, in which a s ° ft appears as a factor. 

It is important to realize that in the |cc~8g) Fock state not only the cc pair is in a colour- 
octet state, but also the three particles, c, c and g, are in an S-state. The covariant spin 
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wave functions we are using for the Xcj read 

Sg> = + M ) (p„/M - 7,), 4? = 7f + M 2) (14) 
The colour-octet component of the Xcj is given by 

lxS> = f /f / dzKka^Czx, z 2 , z 3 ) Sfi, (15) 

where i = A/2 is the Gell-Mann colour matrix and a the colour of the gluon. Since 
orbital angular momenta are not involved, the transverse degrees of freedom are already 
integrated over. Therefore, we only have to operate with a distribution amplitude that 
is, as usual, subject to the condition / dzxdz^j = 1, and an octet decay constant /} . 
In the following, we will make the plausible assumption that the colour-octet Xcj states 
differ only by their spin wave functions, i.e. the distribution amplitudes as well as the 
decay constants are assumed to be the same for the two Xcj states. 

Usually higher Fock state contributions to exclusive reactions are suppressed by powers 
of 1/Q 2 [[TJ where Q = m c in our case. However, for the decays of P-wave charmonia this 
suppression does not appear as a simple dimensional argument reveals: the colour-singlet 
and octet contributions to the decay amplitude defined in @ behave as 

M? ~ flffm-^. (16) 

The singlet decay constant, fj, represents the derivative of a two-particle coordinate 
space wave function at the origin. Hence it is of dimension GeV 2 . The octet decay 
constant, /} , as a three-particle coordinate space wave function at the origin, is also of 
dimension GeV 2 . Since, according to @, Mj is of dimension GeV, n c = 3 in both cases. 
We note that the Xcj decay constants may also dependF] on m c . 

We estimate the colour-octet contributions to the Xcj decays by calculating the set of 
Feynman graphs of which representatives are shown in Fig. |2|. The graphs of group 11 
only contribute to decays into 7r°7r°. In order to simplify matters we carry through the 
calculation in the collinear approximation and assume a 5-function-like ccg distribution 
amplitude (zi = z 2 = (1 — z)/2; z 3 = z ~ (1 — 2m c /M c ) ~ 0.15, where M c is the average 
charmonium mass and z% is the momentum fraction carried by the constituent gluon). 
Most of the details of the calculation, as well as the application of the modified HSA, are 



left to a forthcoming paper p2[ . We note that in some of the graphs shown in Fig. g the 
virtuality of the quark propagator adjacent to the Xcj constituent gluon is small. It is 
these diagrams that constitute, to leading order in a s and z, the higher Fock state |qqg) of 
the pion. However, for arbitrary z, these diagrams alone do not lead to a gauge-invariant 
amplitude. 

In the calculation of the colour-octet contribution, a number of singular integrals 
appear as a consequence of the collinear approximation. If the transverse momenta are 
kept, all integrals are well defined. The difficulty arises from the propagator D ca = 
[(z — x)(z — y) of the gluon exchanged between the light quarks for which the usual 

xe prescription fails. In order to regularize these integrals we replace D ca by 

D= [(z-x)(z-y) + q 2 + te}-\ (17) 



In the XcO case for instance the spin wave for the cc Fock state reads S5 = \p + Mq + 2ff\ If. 

dfcfc 4 1 
167r a Af 



The decay constant is then defined as / (1) = ^ / d ^ ®o\ k ) = -t\R'p{0)\/^/167rm c . 
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T(Xc.J - vr°7r°) [keV] 


BR( Xc j - ttV) [%] 


B 2 = 


25.7 


1.81 


0.18 


0.091 


Exp. PDG [16 


42 ± 18 


2.2 ±0.6 


0.31 ±0.06 


0.110 ±0.028 



Table 3: Decay widths and branching ratios of XcJ ~* 7r°7r° (colour-octet contributions 
included; m c = 1.5 GeV, A QCD = 0.2 GeV). 



The integrals can then be worked out straightforwardly ||23|| . The regulator g represents 



a mean transverse momentum of the quarks inside the pions. Indeed, 

g 2 = (kl)/(4ml) 



We checked that for appropriate values of g this regularisation recipe provides, to a 
reasonable approximation, similar numerical results as one obtains when all transverse 
quark momenta are kept and, weighted by Gaussian /^-dependences, integrated over. 

The final results for the colour-octet contribution to the decay amplitudes can again 
be written in the form ( |I0D with 

a ( j } = mj(z, g) , by = nb, } (z, g) , c { j ] = kcj{z, g) , (19) 

where cij, bj and cj are only weakly dependent on the exact values of the parameters z 
and g. Thus, it appears reasonable to take k, defined by 

« = V^/o'V ( 2 °) 

as the only fit parameter. Evaluating the colour octet contribution from the asymptotic 
7r wave function (B 2 = 0) with z = 0.15, g about 0.1, which corresponds to some typical 
transverse momentum of 300 MeV, and, guided by our results for the singlet contribution 
obtained within the modified HSA, a s (fi e ^ d ) = 0.45, we find 

a = (151.1 + «25.9) GeV~ 2 , a 2 = (156.1 ± z6.6) GeV" 2 . (21) 
The quantities bj and cj, necessary for B 2 ^ case, will be given in |22| . 



The Xcj decay widths into pions are now obtained by adding coherently the respective 
colour-singlet and colour-octet contributions. Since both contributions have the same 
scaling with 1/c they have to be evaluated in the same scheme, i.e. the colour-singlet 
result in the collinear approximation has to be used. The four widths, Xcj — * 7r + 7r~, 
7r°7r° for J = 0, 2 are given in terms of the single non-perturbative factor k. A fit yields 
k = 0.161 GeV 2 and the individual widths are shown in tables and |^. We present also 
predictions for the respective branching ratios using the experimentally measured total 
widths. As shown in tables ^| and [3] the inclusion of the colour-octet mechanism brings 
predictions and data in generally good agreement. 

Finally we address the question whether the value found for the parameter k is sensible. 
For this purpose we estimate the probability for the Xco to be in the colour-octet Fock 
state from a plausible wave function: 

^o\zi, k ±i ) = Nz x z 2 zl exp {-2a 2 m 2 [(z 3 - zf + ( Zl - z 2 f] } 

xexp{-a x 5>L} . (22) 
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This ansatz combines the known asymptotic behaviour of a distribution amplitude for 



a qqg Fock state |24j] with mass-dependent exponentials and a Gaussian fc^-dependence 
in analogy to the Bauer-Stech-Wirbel parametrization of charmed-meson wave functions 
p5| . The mass- dependent exponentials guarantee a pronounced peak of the distribution 



amplitude at Z\ ~ z 2 — m c /M c . The 5-function-like distribution amplitude used in the 
estimate of the colour-octet contribution appears as the peaking approximation to this 
function. Since the ccg Fock state is an 5*-wave state we assume its radius to be equal to 
that of the S'-state charmonia, namely 0.42 fm |15| . In this case the oscillator parameter 



a x takes the value 1.23 GeV . The probability of the colour-octet Fock state is 

Pccg = (/o (8) /2.1 x 10- 3 GeV 2 ) 2 . (23) 



This relation is obtained with z = 0.15; it is however practically independent of the exact 
value of z. Taking a s s oit ~ vr we obtain from (|2§) /,S 8) ~ 0.9 x 10~ 3 GeV 2 and P C c g ~ 19%. 
Hence, the probability for the Xcj to be in the ccg state comes out much more reasonable 
than one might have hoped for. 

In summary, we have shown that the constraints on the pion-distribution amplitude 
from F nl {Q 2 ) data lead to predictions for the exclusive decay Xcj kit that fall well (a 
factor of ten) below the data when restricting to the dominant, colour-singlet |cc( 3 P/)) 
Fock state. The discrepancies persist after the inclusion of transverse degrees of freedom 
and Sudakov suppressions in the modified HSA approach. The deviations are, in fact, so 
large that they cannot sensibly be explained by higher-order perturbative QCD (a s ) or 
genuine higher-twist (l/m c ) corrections. 

We have shown that contributions to the exclusive decay from the 0(w)-suppressed 
Fock component |ccg) are actually of the same power in l/m c as the conventional colour- 
singlet contribution. Predictions including this colour-octet mechanism basically require 
only a single additional non-perturbative parameter. With quite a reasonable value for 
it, we obtain decay widths and branching ratios for the exclusive decays Xcj —* vr + 7r _ and 
7T°7r which are in good agreement with data. 

We do not expect a significant colour-octet contribution to J/-0 (and T) decay into 
pions, since contributions from the higher- Fock state |cc( 3 Pj)g) are truly suppressed by 
1/m 2 . The solution to a similar discrepancy between data and theory for J/tp — > 7r + 7r~ 
is therefore still open and needs further study. Our findings for exclusive Xcj apply 
to P-wave bottomonium decays Xbj ~~ > 7171 as well: a large fraction of the width will 
originate from the bbg state. One may speculate that the colour-octet mechanism is 
relevant to other exclusive decays as well, e.g. B — > XcjK. Work along these lines is in 
progress. 

Acknowledgements: One of us (G.A.S) wants to thank the Theory Group of SLAC for 
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